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CN Abstract 

>• Recently, we introduced a new procedure for computing a class of one-loop BPS-saturated 

K/i amplitudes in String Theory, which expresses them as a sum of one-loop contributions of all 

;_ I perturbative BPS states in a manifestly T-duality invariant fashion. In this paper, we extend 

this procedure to all BPS-saturated amplitudes of the form Jj- r^j^]^^^ 0, with €> being a weak 
(almost) holomorphic modular form of weight —k/2. We use the fact that any such €> can be 
expressed as a linear combination of certain absolutely convergent Poincare series, against 
w^hich the fundamental domain T can be unfolded. The resulting BPS-state sum neatly 
exhibits the singularities of the amplitude at points of gauge symmetry enhancement, in a 
chamber-independent fashion. We illustrate our method with concrete examples of interest 
in heterotic string compactifications. 
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1 Introduction 

Scattering amplitudes in closed string theory involve, at ?2-th order in perturbation theory, an 
integral over the moduli space of conformal structures on genus h closed Riemann surfaces. 
The torus amplitude (corresponding to /z = 1) is particularly relevant, as it encodes the 
perturbative spectrum of excitations. Moreover, for special choices of vacua and of external 
states, corresponding to a special class of f -term interactions in the low energy effective 
action, the torus contribution exhausts the perturbative series, and thus can serve as a basis 
for quantitative tests of string dualities (see e.g. [1 [ and references therein). 

The moduli space of conformal metrics on the torus is the Poincare upper half plane H, 
parameterised by the complex structure parameter t = Ti + 1x2, modulo the action of the 
modular group SL(2, Z). After performing the path integral over the world-sheet fields and 



over the location of the vertex-operator insertions, the relevant amplitude is then expressed 
as a modular integral 

j d^A{Tl,T2), (1.1) 

where J-" = {r € H | — ^ < ti < j, |t| > 1} is the standard fundamental domain, d}i = 
T2"^ dzi dT2 is the SL(2, Z)-invariant integration measure, and ^ is a modular-invariant func- 
tion whose precise expression depends on the problem at hand. With this choice of integra- 
tion domain, the imaginary part T2 can be identified with Schwinger's proper time, while 
the real part Ti is the Lagrange multiplier imposing the level-matching condition. Part of 
the difficulty in evaluating integrals of the form ( |1.1| | is the unwieldy shape of T, which 
intertwines the integrals over Ti and T2. 

Depending on the function A{ti, T2) methods have been devised to overcome this prob- 
lem. If ^ is a weak almost holomorphic functioiylof r (or, alternatively, an anti-holomorphic 
function), the surface integral over J-" can be reduced by Stokes' theorem to a line-integral 
over its boundary dT that can be explicitly computed [2|. On the contrary, if ^ is a gen- 
uine non-holomorphic function, as is the case for the one-loop partition function of closed- 
oriented strings, no useful method is know^n to evaluate the integral, but one can use the 
Rankin-Selberg-Zagier transform \3] to connect the integral to the graded sum of physical 
degrees of freedom |ffll7|. A frequently encountered intermediate case is that of modular 
integrals of the form 

I d^^a+kAG,B,Y;T^,T2)0{T), (1.2) 

w^here 

Fd+kAG, B, Y; Ti, T2) ^ rf 2 ^ c^lvl -^\vl (I.3) 

Pl,Pr 

is the partition function of the Narain lattice of Lorentzian signature {d + k,d), G, B, Y pa- 
rameterise the Narain moduli space SO(d + k,d) /SO{d + A;) x SO(d), and <E'(t) is a w^eak 
almost holomorphic modular form of negative w^eight w = —k/2, which we shall refer to as 
the elliptic genus. Such integrals occur in particular in one-loop corrections to certain BPS- 
saturated couplings in the low energy effective action of heterotic or type II superstrings. 

The traditional approach in the physics literature for computing modular integrals of the 
form ( |1.2| has been the orbit method, which proceeds by unfolding the integration domain 
J-" against the lattice partition function T^j+^d ||8-19|. While this procedure yields an infi- 



nite series expansion which is useful in certain limits in Narain moduli space, it does not 
make manifest the invariance under the T-duality group 0(d + A:, d, Z) of the Narain lattice, 
nor does it clearly display the singularities of the amplitude at points of gauge symmetry 



^ By iveak almost holomorphic we mean an element in the graded polynomial ring generated by the holomorphic 
Eisenstein series £4 and Eg, the almost holomorphic Eisenstein series £2 and the inverse of the discriminant 1/ A. 
Our notations for Eisenstein series and other modular forms are collected in Appendix |A.1| The adverb weak 
refers to the fact that the only singularity is, at most, a finite order pole at the cusp i^ = 0. 



enhancement. 

In 1 20 1 we proposed a new method for dealing with modular integrals of the form \i.2) , 
which relies on representing the elliptic genus as a Poincare series, and on unfolding the 
integration domain against it rather than against the lattice partition function. The advan- 
tage of this method is that T-duality remains manifest at all steps, and the result is valid 
in all chambers in Narain moduli space, unlike the conventional approachrl Moreover, the 
amplitude is expressed as a sum over all BPS states in the spectrum, thus generalising the 
constrained Eisenstein series constructed in [|2 1 J. p] Finally, the singularities of the amplitude 
at points of enhanced gauge symmetry can be immediately read-off from the contributions 
of those BPS states w^hich become massless. 

The main difficulty in implementing this strategy is due to the fact that the standard 



Poincare series representation of a w^eak holomorphic modular form of weight w < |28- 
Soj is only conditionally convergent, and therefore unsuited for unfolding. In f20l we at- 
tempted to circumvent this problem by considering a class of non-holomorphic Poincare se- 
ries £(s, K, w) that provide a natural regularisation of the modular forms of interest by insert- 
ing a Kronecker-type convergence factor t^^'" in the standard sum over images. Therefore, 
the resulting Poincare series, originally studied in [31 1, converges absolutely for 3?(s) > 1, 
and the modular integral J^r^jrf£(s,K, w) can be computed by unfolding J^ against it, at 
least for large s. The result should then be analytically continued to the desired value s = j, 
w^here E(s, k, w) becomes formally a holomorphic function of t. This procedure would then 
allow to compute the modular integral ( |1.2| | for any (P which can be expressed as a linear 
combination of such £(^,K,a;)'s, at least in principle. However, this strategy turned out to 
be quite difficult in practice, since this analytic continuation depends on the notoriously sub- 
tle analytic properties of the Kloosterman-Selberg zeta function which appears in the Fourier 
expansion of E(s, k, w). That is the reason why the analysis |20| was restricted to the case of 
zero modular weight, where the analytic continuation is fully under control. 

In the present w^ork, we overcome these difficulties by employing a different class of non- 
analytic Poincare series introduced in the mathematics literature by Niebur |[32| and Hejhal 



| [33| and studied more recently by Bruinier, Ono and Bringmann |34]-37j. Similarly to the 
Selberg-Poincare series E{s,k,w), the Niebur-Poincare series J^{s,k,w) converges absolutely 
for !R(s) > 1, and formally becomes holomorphic in t at the point s = j. However, the 
Niebur-Poincare series can be specialised to the other interesting value s = 1 — j, which 
lies inside the domain of absolute convergence when the weight w is negative. Although at 
this value ^{s, k, w) belongs to the more general class of weak harmonic Maass forma^ that 
are typically non-holomorphic functions of T, it has the important property that any linear 
combination of J-'(l — ^, k, w), whose coefficients are determined by the principal part of a 
given weak holomorphic modular form <J>, is actually a weak holomorphic modular form. 



^See for instance |18| for a detailed discussion on chamber dependence of the traditional unfolding method. 

^BPS states sums have appeared in earlier works 1 22 -26 J. In our approach these BPS sums follow directly 
from unfoldir\g the fundamental domain against the elliptic genus, without any further assumption. 

*A harmonic Maass form is an eigenmode of the weight-w Laplacian on H with the same eigenvalue as weak 
holomorphic modular forms. The positive frequency part of a weak harmonic Maass form is sometimes known 



as a Mock modular form. See Section 2.3 for a more precise definition of weak harmonic Maass forms. 



and equals <P itself. Therefore, given any weak holomorphic modular form 0, the integral ( |1.2| | 
can be computed by decomposing into a sum of Niebur-Poincare series, and by unfolding each of 
them against the integration domain. Moreover, the same strategy works also for weak almost 
holomorphic modular forms (i.e. involving powers of E2)/ where now one has to specialise 
the Niebur-Poincare series to the values s = 1 — ^ + n, w^ith n a non-negative integer. 

The outline of this work is as follow^s. In Section|2| we introduce the Niebur-Poincare se- 
ries J^{s,K,w), discuss their main properties, present their Fourier coefficients and identify 
their limiting values at s = 1 — ^ + n. We conclude the section by showing the important 
result that any weak almost holomorphic modular form can be represented as a linear combi- 
nation of them. In Section 3 we evaluate the modular integral fj- F^^^^i^ T{s, k, — j ) in terms 
of certain BPS-state sums and discuss their singularity structure. In Section [4| we use this 
result to compute a sample of modular integrals of physical interest of the form ( |1.2| |. In 
Appendix [A] we define our notation for modular forms, we collect various definitions and 
properties of Whittaker and hypergeometric functions, and we introduce the Kloosterman 
sums and the Kloosterman-Selberg zeta function. Finally, in Appendix Ib] we briefly discuss 
the relation between the Selberg- and Niebur-Poincare series, and betw^een the "shifted con- 
strained" Epstein zeta series and the above BPS-state sums. The reader interested only in 
physics applications may skip Section|2]and proceed directly to SectionlSl which begins with 
an executive summary of the main properties of J^{s, k, w). 

Note. After having obtained most of the results in this paper, we became aware of ref. 
p4| where similar computations have been performed for general even lattices of signature 
(d + k,d) with d = 0,1,2, in particular reproducing Borcherds' automorphic products for 
d = 2 f38|. Unlike fsl^, we restrict the analysis to even self-dual lattices (with k = mod 8), 
which suffices for our physics applications, but we allow for almost holomorphic modular 
forms and arbitrary dimension d. 

2 Niebur-Poincare series and almost holomorphic modular forms 

In this section, we introduce the Niebur-Poincare series T{s, k, w), a modular invariant reg- 
ularisation of the naive Poincare series of negative w^eight. We present its Fourier expansion 
for general values of s, and analyse its limit ass— >1 — ^ + n where n is any non-negative in- 
teger. We explain how^ to represent any w^eak almost holomorphic modular form of negative 
weight as a suitable linear combinations of such Poincare series. 

2.1 Various Poincare series 

In order to motivate the construction of the Niebur-Poincare series, let us start with a brief 
overview of Poincare series in general. Let w be an even integeijj and / a function on the 
Poincare upper half plane H. The action of an element 7 = ( " , ) G T = SL(2, Z) on / is 

^In this paper we shall restrict to the case of even weight w in order to avoid complications with non-trivial 
multiplier systems, though the construction can be generalised to half-integer weights. 



given by the Petersson slash operator 

(/U7)(T) = (cT + d)-/(7-T), 7.T=^. (2.1) 

If / is invariant under Too = ( „ *A C F, the Poincare series of seed f and weight w 

P{f,w;T)^P{f,w) = l Yl fl-^-y (2-2) 

defines an automorphic form of weight w on H, which is absolutely convergent provided 

1— - 
/(t) <C Tj ^ as T2 — ?► 0. As an example, the choice /(t) = q^''^ with w > 2 leads to the usual 

holomorphic Poincare series 

P{k,zv) = 1 ^ (cT + d)-"'e-2^^'^iT3, (2.3) 

{c,d)=l 

w^here the pair {a,b) is determined modulo {c,d) by the condition ad — be = 1. Depending 
on the value of k, Eq. ( |2.3| describes different types of modular forms. For k = 0, P{k, w) is 
actually an Eisenstein series, while for k < —1 it is a cusp form, and must therefore vanish 
if 2 < w < 12, an observation that will be important later. For k > 0, Eq. ( |2.3[ ) represents 
instead a weak holomorphic modular form w^ith a pole of order k at q = 0, P{k, w) = q^^ + 
0{q). 

For w < 2, the Poincare series ( |2.3| l is divergent and thus needs to be regularised. One 
possible regularisation scheme, introduced in the mathematical literature in |28 29| and dis- 



cussed in the physics literature in |30 1, is to consider the convergent sum 



P(k,w) = \ Wm y y (cT + d)-^V2^i'^?SlR( ^^^^) , (2.4) 

'^--Ic^K \,\<kc4U UcT + d)) 

where i? is a specific regulating factor such that R{x) ~ x^~^" /r{2 — a;) as x — )■ and 
approaches 1 as x — ?► oo. While this regularisation preserves holomorphicity, it does not nec- 
essarily produce a modular forrrr] except for small \w\ where the modular anomaly can be 



shown to vanish. Moreover, the convergence of (2.4 1 is conditional, which makes it unsuit- 
able for the unfolding procedure. 

Another option, introduced by Selberg pTj and considered in our previous work pO| , is 
to jettison holomorphicity and introduce a convergence factor a la Kronecker, thus consider- 
ing the Poincare-series 

£(s,K,«;)^l E r-^^^{cT + d)-e-'"^^l^^ (2.5) 

^The holomorphic Poincare series \2A) is in general an Eichler integral, i.e. a function F(t) which satisfies 

F(t) — (F|j(,7)(t) = rj{T) where r^ is a polynomial of degree —iv in t, whose coefficients depend on a,b,c,d. 



We shall comment in Section 2.3 on the modular completion of P{k, w) 



s— 



associated to the seed /(t) = Tj ^ q ^. We shall refer to ( |2.5| as the Selberg-Poincare series. 



The series (2.5 1 converges absolutely for K(s) > 1 and becomes formally holomorphic at 
s = ^■. However, for w < 2 this value lies outside the convergence domain, and the ana- 
lytic continuation to s = f depends on the analytic properties of the Kloosterman-Selberg 
Zeta function, defined in Appendix |B| which are notoriously subtle. In particular this an- 



alytic continuation generally leads to holomorphic anomalies. For this reason, in |20| w^e 
restricted the analysis to the case a; = 0, where the analytic continuation is under control. 
Another drawback of the Selberg-Poincare series ( |2.5| l is that it fails to be an eigenmode of 
the Laplacian on H, rather it satisfies [391 

[^lo + 2s(l-s) + \w{w + 2)\ E{s,K,w) =271k{s- f)E(s + l,?c,ii;), (2.6) 

where A^i, is the weight-ro hyperbolic Laplacian defined in ( |A.1| |. Since £(s + 1, k, w) may in 
general have a pole at s = |^, the analytic continuation of £(s, k,w) to this value is not even 
guaranteed to be harmonic. 



To circumvent these problems, following 1 32 - 34 1 we introduce a different regularisation 



of the Poincare series ( |2.3| l for negative w^eight, w^hich is both modular invariant and an- 
nihilated by the operator on the l.h.s. of ( |2.6[ ). Namely, we choose the seed in ( |2.2| to be 

/(t) = Ms,zoi-KT2) e-2^i'^'^i where 

MsAy) = |47ryr? M.3g,(y)^,_i (47r|y|) (2.7) 

is expressed in terms of the Whittaker functioirlMA,^ (z) . We thus define the Niebur-Poincare 
series 



7G-r„\r 

KT2 \ f ^, fa CTi+d 

C cIcT + dp 






{c,d)=l 



(2.8) 



Since A4s,w{y) ~ i/^(^)^l as y — > 0, Eq. ( |2.8| converges absolutely for K(s) > 1, indepen- 
dently of w and k. Moreover, for k > 0, the case of main interest in this work, the seed 
behaves as 

MsA-^T2)e-^-''-^ ^ -^^^q-' as T2 ^ oo , (2.9) 

i [s + 2) 

so that T{s,K,w) can indeed be viewed as a regulated version of the naive Poincare series 
P{q^'^, w), up to an overall normalisation. By construction it is an eigenmode of the weight-if 
Laplacian on H, 

[Arv + ls{l-s) + lw{zv + 2)] F[s,K,w) =0, (2.10) 



For a definition of Whittaker functions and some of their properties see Appendix 



A.2 



for all values of s, k, w. We shall denote by 7^(5, w) =7^(1 — 5, w) the space of real-analytic 
solutions to ( 2.10| | which transform with modular weight w under F. 



The raising and lowering operators D^y, Dju defined in ( A.2| |, map T-L{s, w) into T-L{s,w i: 
2), and have a simple action on the Niebur-Poincare series 



D-,„-J^{s,K,w) =2k{s + '^)T{s,k,w + 2), 

1 ,„_. _. (2-11) 

8k 



D-,„-J^{s,K,zv) = — (s- f ) J'(s,K,a;-2) 



Furthermore, under the action of the Hecke operator ( A.7| J-^(s, k, w) transforms as 



T^, ■ T{s,K,w) = J^ d^-'"J^{s,KK'/d^,w). (2.12) 

d\(K,K') 

In particular, setting k = 1, the series J^{s, k', w) is obtained by acting with T^' on J^{s, 1, w ). 



While the Poincare series (2.8 1 converges absolutely only for K(s) > 1, it is know^n to have 
a meromorphic continuation to the complex s-plane, holomorphic in the region !K(s) > j 



1 32 40 , but with poles on the lines s S ^ + iR and s G ^ + iR. Moreover, the 'completed' 



series 



is known to be invariant under s i— > 1 — s, up to an additive contribution proportional to 



the non-holomorphic Eisenstein series E*{s,zv) 1 32 401. In this work however we shall only 
consider T{s, k, w) in its domain of convergence !K(s) > 1, except for w = where we allow 
s = l. 

2.2 Fourier expansion of the Niebur-Poincare series 

The Fourier expansion of J^{s,k,w) can be obtained follow^ing the standard procedure of 
extracting the contribution from c = 0,d = 1, setting d = d' + mc in the remaining sum, and 
Poisson resumming over m. The result is p4}[36) 



T{s,K,w) = MsA-KTi) e-^""'"^' + Y^ ^m{s,K,w) e^""'"'^^ , (2.14) 



w^here, for zero frequency 



-s-^ 



^, . 22— i— 7ri+«-2K«-2r(2s-l)c7i_2sW 1-.-^ ,_^, 

^'^'' ""' ^^ = r(s-f)r(s + f)a2s) ^- ' ^^-''^ 



w^hile for non-vanishing integer frequencieaj 

47rKi"'T(2s) \m 
r{s + f sgn(m)) I K 



47rKi"'T(2s) imif 
J'm{s,K,w) = ———^^ T^ ~ 2^s{m,-K)Ws,iuimT2). (2.16) 



'^Note that J^-k<o does not include the contribution from the first term in (2J4b 



In these expressions, (Xg (/c) = J^m d^ is the divisor function and Zs {in, —k) is the Kloosterman- 
Selberg zeta function ( |A.37[ ), a number-theoretical function which plays a central role in 
the theory of Poincare series. The function yVs,w is expressed in terms of the Whittaker W- 
function as 

WsAy) = \^rzy\-'i W.,g^(^)^,_, {An\y\) , (2.17) 

and is determined uniquely by the requirement that Ws,zu(nT2) e^/rimTi ^^^ annihilated by the 
Laplace operator on the l.h.s. of ( |2.10| |, and be exponentially suppressed as T2 ^ oo. 

Using the properties ( A.35| l and ( A.36| l, it is straightforward to check that all Fourier 



modes transform according to (2.11 1 under the raising and lowering operators Dw, Du,. More- 
over, using the action ( |A.9| | of the Hecke operators on the Fourier coefficients, and the Selberg 
identity ( |A.39| | satisfied by the Kloosterman sums, one can show that 

T„ ■ J'(s, 1, w) = J'(s, K, w) . (2.18) 

Eq. ( 2.12| follows then from this equation and from the Hecke algebra ( A.8 1. 



2.3 Harmonic Maass forms from Niebur-Poincare series 

Let us focus on the Niebur-Poincare series J^{s,k,w) at the point s = 1 — ^. To motivate 
this value, we recall that any weak holomorphic modular form is an eigenmode of Aa; with 
eigenvalue — f , and therefore belongs to H{s,w) for s = 1 — ^ (or equivalently, s = j). 
However, w^eak holomorphic modular forms are not the only eigenmodes of A^ with this 
eigenvalue. In fact, the space T-L{1 — ~^,w) is known as the space of weak harmonic Maass 
forms of weight w, of which weak holomorphic modular forms are only a proper subspace. 

The Fourier expansion of a general w^eak harmonic Maass form of w^eight w is given 
by EI) 



))(=— 00 



1 ^- ^4:71'^-^^^'^ 

W 



^= E (-m)^"-! ^_,„ r(l - w, -AKmT2) c,'^ + ^° ^^""^l + E «-^'"' (2-19) 



where r(s, x) is the incomplete Gamma function and am,bm are coefficients constrained 
by modular invariance. As a result, a generic weak harmonic Maass form has an infi- 
nite number of negative frequency components, which are non-holomorphic functions of 
T. A harmonic Maass form splits into the sum <J> = cpQ + <J>(, of a holomorphic part 0a = 
'}2m=-K '^m ^"\ sometimes called a Mock modular form, and a non-holomorphic part 0j,. The 
non-holomorphic and holomorphic parts can be extracted using the lowering operator D^, 
and the iterated raising operator D^^^". Indeed, 

• the operator D^, annihilates the holomorphic part, and produces, up to powers of T2, 
the complex conjugate of a holomorphic modular form Y of weight 2 — w, 

D,,-0 = A„ ■ 0, = -2I-2- (7rT2)2— W , Y{t) = E b,„ q"" , (2.20) 

m=0 



sometimes known as the shadow. 

• the iterated raising operator Dl^^, also known in the physics literature as the Farey 
transform \42\ , annihilates the non-holomorphic part, and produces a weak holomor- 
phic modular form S of weight 2 — w, 

oo 

Di-.cD = Di,--0„ = S, S^ E (-2m)i-fl„^-, (2.21) 

that we shall call the ghost. The ghost encodes the holomorphic part of the harmonic 
Maass form (modulo an additive constantjj 

Returning to the Niebur-Poincare series, we see that by construction the series J^{s, k, w) 
at the special point s = 1 — j — which, for iv < 0, belongs to the convergence domain — is 
a w^eak harmonic Maass form of weight w. Indeed, using ( |A.30 | we find that its Fourier 



expansion (2.141 reduces to 



J-(l - f ,K,z(;) = A1i_|,,„(-KT2) e-2--i + Y, -^m(l - f ,K,a;) e^-'"-! , (2.22) 

meZ 

w^here the seed simplifies to a finite sum 

-Mi-f ,.(-KT2) e-2--i = r(2 - w) L-^ - r f (4^^^ ^^ ^^^ 

= [r{l-w) + {l-w)r{l-wAnKT2)]q~'', 
and the remaining Fourier coefficients reduce to 



-m '^ 2 

K 

r2. 



Tm<o{l-^_,K,w) =Anr'"K{l-w)i^^\ Z^_m{m, -k) r{l - w, -AnmT2) e-^'""^^ 



(2.24) 

^Notice that the ghost is only defined for integer weight w, unlike the shadow, which extends to the case of 
half-integer weight Mock theta series. 



One thus recognises an expansion of the form ( 2.19| l with coefficients 



fl_K = r{2-zv), 



flo 



flm>o=47ri ^"Kr{2 — w) ( — ) Z^w{m,—K), 






(2.25) 



bo = 0, 

&„,>0 = (1 -w)K^''"SniK+4^ni^" {1 -w){mK)^ 2 Z^_a;(m,K). 

2 

In particular, &o = 0/ so that the shadow of J-'(l — ^,K:,zf ) is a cusp form of weight 2 — w, 
proportional to the holomorphic Poincare series P{—k,2 — w). Indeed, using ( 2.11| we find 



1 — m 

D,„ ■ T{1 - l,K,w) =^-T{\ - l,K,w- 2) 



-, _ (2.26) 

^(47rKT2)2-"'P(-K,2-a;), 



8k 

where in the second line we have recognised the Fourier expansion of the standard holo- 
morphic Poincare series of weight greater than 2. Similarly, using ( 2.11| the ghost of J-'(l — 
^^,■K,w) is 

Bl-'" ■ T(\ - f , K, w) = (2k)1-'" T(2 - w) T{\ - f , k, 2 - w) , (2.27) 

and corresponds to the Niebur-Poincare series J-'(^,K,if'), with w' = 2 — w > 2 within 
the convergence domain. Moreover, the Fourier expansion of the latter reproduces that of the 
Poincare series V{%,w^ of positive weight 

^(^,K,z.')=r^ + 27ri-' E (-)^E^^^^^^^^wf-v/-.)^'". (2.28) 

As an aside, we note that the holomorphic part ^-'^(l — ^,k,w) of the Niebur-Poincare se- 
ries J^{s,K,w) at s = 1 — ^ reproduces the Fourier expansion of the Poincare series r(2 — 



iv)P{k,w) defined by holomorphic regularisation as in (2.4| and worked out in 1 28 291. 
Therefore, the non-holomorphic part Ji,(l — j,k,w) of the same Niebur-Poincare series 
provides the modular completion of the Eichler integral P{k,w) — a clear advantage of 
modular-invariant regularisation over holomorphic regularisation. 

To make this discussion less abstract, we shall now exhibit the harmonic Maass form 
^(1 — j,K,w), its shadow and its ghost for the two cases w = —10 and w = —14 (lower 
values of |w| will be discussed in the next subsection) and k = 1. Evaluating the Fourier 
coefficients numerically, we find: 



10 



For w = —10, 

J-(6, 1, -10) = J^b{6, 1, -10) + 11! 



,-1 



65520 
691 



1842.89^-23274.08^^ 



(2.29) 



■where Ji,(6, 1, —10) is the non-holomorphic component. The shadow of ( 2.29| l reads 

J^{6, 1, -12) = (47rT2)i2 P(-l,12) , P{-1, 12) = ^12 ^ , 



(2.30) 



where the modular discriminant A generates the space of cusp forms of weight 12. The 
ghost, obtained by acting with D^^ on the holomorphic part, can be written as 



^(6,1,12) 



83£3E? - 11£| 



■^4'^6 



+ OC12 



(El - Elf 
72A ' "'^ A 

q-^ + 1842.89 q + 47665306.53 q^ 



(2.31) 



where ^12 = 0.201029508104 ...,,612 = 2.840287517 ... are irrational numbers. The co- 
efficients 1842.89, 23274.08, 47665306.53 are two-digit approximations of the exact val- 
ues 324(9216*12 - 1847), (60617- 69984 ai2)/2, 1024(60617 -69984 0:12), respectively. 



This example was discussed in detail in |37|. 
Similarly, for w = — 14, 



J-(8, 1,-14) =J-fo(8, 1,-14) + 15! 



1 16320 9 

r - ^^ - 45.67^ - 366.47^2 ^ 



(2.32) 



w^here ^j,(8, 1, —14) is the non-holomorphic component. The shadow of ( 2.32| | reads 
J-(8, 1,-16) = (47rT2)i*^P(-l,16), P(-l,16) = ^leE^A 



(2.33) 



where E^A generate the space of cusp forms of weight 16. The ghost, obtained by acting 
with D^^ on the holomorphic part, can be written as 



7oc:4p2 r p4 p7 of4i72 i^ p p4 

J- (8, 1, 16) = ^^^ h ocie ^ 

= q-^ + 45.67 q + 12008361.57 (j^ + . . . , 



(2.34) 



where ai6 = 0.137975847804... and j6i6 = 1.3061364711... are irrational numbers. 
The coefficients 45.67, 366.47 in Eq. I\2.32) are two-digit approximations of the exact 
values 36(82944ai6 - 11443), (314928 aie - 37589)/16, respectively 

These two examples illustrate the fact that Fourier coefficients of harmonic Maass forms 
are in general irrational numbers. 
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w 


-^(l-f,l,^i^) 


T{l-l,l,2-w) 





; + 24 


ElEeA-^ 


-2 


3! £4^6^"^ 


£4(;-240) 


-4 


5! £4^-1 


£6(;+204) 


-6 


7\E(,A-^ 


£|(;-480) 


-8 


9\E^A-^ 


£4£6(; + 264) 


-12 


13\A-^ 


E|£6(; + 24) 



Table 2.1: Weak holomorphic modular forms obtained as the limit s — )■ 1 — ^ of jF{s, \,w), 
for the values w G {0, —2, —4, —6, —8, —12}. For w negative and outside this range, the limit 
yields a weak harmonic Maass form. The second line shows the ghost, which is a weak 
holomorphic modular form of weight 2 — w w^ith vanishing constant term (aside from the 
case w = 0) 

2.4 Weak holomorphic modular forms from Niebur-Poincare series 

We now come to our main goal, i.e. to find an absolutely convergent Poincare series repre- 
sentation of any weak holomorphic modular form 0r„ of weight w < and K-order pole at 
the cusp, with given principal part 

^7o{t)= E «-'?'"• (2-35) 

-K<m<Q 

As we shall see, any such 4>;y can be expressed as a linear combination of the Niebur-Poincare 
series F{s,k,w). 

We have noted in the previous subsection that the eigenvalue of a w^eak holomorphic 
modular form under the hyperbolic Laplacian Afy coincides with the eigenvalue of the Niebur- 
Poincare series whenever s = \ — ^■. At this value, however, J-'(l — ^,k:, zf)isa weak har- 
monic Maass form, in general not holomorphic. Exceptions to this statement occur at the 
special values w € {—2, —4, —6, —8, —12}, where the space of holomorphic cusp forms of 
weight 2 — w is empty, and T{\ — ^^,k,w) can be recognised as an element of the ring of 
weak holomorphic modular forms by matching the principal part of their expansions. For 



K = 1 the exact identification is reported in Table 2.1 while for k > 1, the proper identifica- 
tion of J^(l — J, K, w) can be obtained by acting on J^{1 — j, 1, w) with the Hecke operator 
Tk, as given by Eq. ( |2.18[ >. 

For w < outside the list above, the space of cusp forms of weight 2 — w is not empty, 
and T{1 — j,k,w) is indeed a genuine harmonic Maass form, with non-vanishing shadow. 
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Nevertheless, it can be shown 1 34 1 that the linear combination 

l 

G{s,w) = — r J^ a„iJ^{s,m,w) , 

^ l^ ^) -K<m<0 

with coefficients fl^ determined by the principal part 

-K<m<0 



(2.36) 



(2.37) 



of any weak holomorphic form 0^^ of negative weight w, reduces to a weak holomorphic 
modular form for s = 1 — j, namely (Piu itself. Said differently the shadows of the weak 



harmonic Maass forms J-'{1 — j,k,zv) cancel in the linear combination (2.361. As a result, 
any 0j„ can be represented as the linear combination 



'Pm 



^ I —K<m<0 



or equivalently using (2.231, as an absolutely convergent Poincare sum 



0u 



7Gr„\r 



, -K<m<Oe=0 



1, 



(2.38) 



(2.39) 



where the subtraction in the bracket ensures that the seed is ©(Tj ^ ) as T2 — > 0. We stress 
that, unlike the holomorphic regularisation in ( |2.4| |, the expression ( |2.39| l is manifestly mod- 
ular covariant and absolutely convergent. 

To illustrate the power of Eq. ( 2.38| l, let us reconsider the two examples of the previous 
subsection, now^ allow^ing for k = 1,2. 

• For w = —10, J-'(6,2, — 10) and ^(6,1,-10) are separately weak harmonic Maass 
forms with irrational coefficients, but the linear combination 



hf, hi 



J^(6,2, -10) +24J'(6, 1,-10) = 11!^ = 11! (^"^ + 24 ^"^ - 196560 + ... ) (2.40) 

produces (up to an overall normalisation) the unique weak holomorphic fornv^ of 
weight —10 with a double pole at ^ = 0; 

Similarly, for w = —14, J^(8,2, —14) and T{S, 1, —14) are separately weak harmonic 
Maass forms with irrational coefficients, but the linear combination 



J"(8, 2, -14) - 216 J'(8, 1, -14) = 15! ^ = 15! {q-^ - 216 q-^ - 146880 



(2.41) 



produces (up to an overall normalisation) the unique w^eak holomorphic form of weight 
— 14 with a double pole aiq = Q. 



10, 



Compare the simplicity of our expression to the corresponding equation in Sec 4.1 of J37J. 
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Similar relations occur for higher negative weight w < —lA and higher order k of the 
pole atq = 0. 

2.5 Weak almost holomorphic modular forms from Niebur-Poincare series 

For physics applications it is important to extend our previous analysis to the case of weak al- 
most holomorphic modular forms, i.e. elements of the ring generated by the almost holomor- 
phic Eisenstein series £2 arid the ordinary weak holomorphic modular forms, or equiva- 
lently, by the modular derivatives D"(P of ordinary weak holomorphic modular forms. 

To this end, it is important to note that for any integer n > 0, it follows from ( |2.11[ ) that 
the Niebur-Poincare series J-'(s, k, w) evaluated at the point s = 1 — ^ + n can be expressed 
as 



1 

(2k)" n! 



^■(1- f + n,K,a;) = ^^^,^„ ^^, D" J"(l - f + n, K,a; - 2n) , (2.42) 



w^here D" is the iterated modular derivative ( A.6| |. The Niebur-Poincare series T{s' ,k,w') 



appearing on the r.h.s. satisfies s' = 1 — ^, and thus is a harmonic Maass form. 

As a result, provided that the coefficients fl„, in the linear combination ( 2.36| l are chosen 
such that 

is the principal part of a weak holomorphic modular form (Pw-in of weight w — 2n, then 
the linear combination Q{s,w) in ( |2.36| l evaluated at the point s = 1 — ^ + n reproduces an 
almost holomorphic modular form of weight w, 

G{l--^ + n,w)= \ Y^ anrJ^{l-'j + n,m,w) = D"0^^2n. (2.44) 

^ V^ ^) -K<m<0 

More generally, we refer to the space 0„>o'H(l — j +n,w) as the space of "weak almost 
harmonic Maass forms", of which almost holomorphic modular forms are only a subspace. 
The general Fourier expansion of such forms can be obtained by taking the limit s = 1 — 
f + n in eqs. ( |2.14[ ) and ( |2.16[ ), and by using the identities ( |A.33| and ( |A.34| |. Similarly, the 



series J^{s, k, w) at the point s = 1 — ^ + n forn < — 2 may be obtained from J-'(s', k, w') at 
the point s' = — ^ by using the lowering operator D-o). 

2.6 Summary 

To summarise this discussion, it is useful to consider the plane of the variables (f ,s) as in 



Figure 2.1 The Niebur-Poincare series J^{s, k, zv) converges absolutely for s > 1. For integer 
values of s, it is generally a weak almost harmonic Maass form, and on the line s = 1 — ^ 
(and w < 0), J^{s, k, w) becomes a weak harmonic Maass form. On the line s = —j, obtained 
from the former by acting with the lowering operator D, T{s, k, w) reduces, up to an overall 
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^ 



7 



ii'. 






■■■x, ^. 



weak almo$t harmonic 
3 



D 

>, " 2 






1 ^3^ 






D / // 

«- • ^ • 

/ :v 



24, ^iEJA 



w 



Figure 2.1: Phase diagram for the Niebur-Poincare series T{s,k,w) for integer values of 
{j,s) with s > 1. For low negative values of w, jF{s,k,w) reduces to an ordinary weak 
almost holomorphic Maass form, see Table 2.2 
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multiplicative factor Tj"'", to the complex conjugate of a cusp form of weight 2 — w, known 
as the shadow of the harmonic Maass form F{s,k,w + 2). On the line s = f , F{s,k,w) 
is instead a weak holomorphic modular form. It is connected to its expression on the line 
s = 1 — ^ by the action of the iterated raising operator D^^'", and thus we refer to it as the 
'ghost' of the harmonic Maass form F{s, k,2 — zv). In the quadrant w > 2,s > 1, F{s, k, w) is 
more generally a weak almost holomorphic modular form. For low negative values of w and 
s integer, J^{s, k, w) is in fact always a weak almost holomorphic modular form, as displayed 
in Table [Z2[ Genuine harmonic Maass forms start appearing at s = 6 and s > 8. 



3 A new road to one-loop modular integrals 

We are interested in the evaluation of one-loop modular integrals of the form ( |1.2| |, while 
keeping manifest at all steps the automorphisms of the Narain lattice, i.e. T-duality. Such in- 
tegrals encode, for instance, threshold corrections to the running of gauge and gravitational 
couplings. The function 0, related to the elliptic genus and dependent on the vacuum un- 
der consideration, is in general a weak almost holomorphic modular form of non-positive 
weight. For example, in J\f = A compactifications of the SO (32) heterotic string (w^ith van- 
ishing Wilson lines) one finds a linear combination of zero-weight w^eak almost holomorphic 
modular forms ||2| 

0{r) = tstrF^ + 2^25 I ^^'^^' + A^^ ''^''^"^' 




E2E^_EM\ ^^^^p2^^j^2 



(3.1) 




2^^^-2'3A ts{trF^)\ 



w^here tg is the familiar tensor appearing in four-point amplitudes of the heterotic string, and 
F and R are the gauge field strength and curvature tw^o-form. A similar expression arises for 
gauge and gravitational couplings in the Eg x Eg heterotic string. 

While the traditional procedure for evaluating integrals of the form ( |1.2| | has been to un- 
fold the integration domain J-" against the lattice partition function T^+if rf, in [|20J w^e instead 
proposed to represent as a Poincare series of the form ( |2.2| , which is then amenable to the 
unfolding procedure. The advantage of this approach is that T-duality is kept manifest at all 
steps and the final result is expressed as a sum over BPS states which is manifestly invari- 
ant under 0{d + k,d;'Z.). Moreover, singularities associated to states becoming massless at 
special points in the Narain moduli space are easily read off from this representation. 

3.1 Niebur-Poincare series in a nutshell 

In order to implement this strategy, it is essential to represent as an absolutely convergent 
Poincare series, so that the unfolding of the fundamental domain is justified. Fortunately, as 
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discussed in detail in Section[2]and summarised in the following, any weak almost holomor- 
phic modular form 0^ of weight w < can be written as a linear combination of Niebur- 
Poincare series, defined as 

Here A4s,iv is related to the Whittaker M-function via 

MsA-y) = (47ry)"'''/'M_.^,_i(47ry), (3.3) 

and s is a complex parameter, the real part of which must be larger than 1 for absolute 
convergence. The choice of the Whittaker function in ( |3.2| | is dictated by the requirement 
that J-'(s, K, w) be an eigenmode of the hyperbolic Laplacian Aa, (see Eq. ( 2.10| ), and behave 



as q^'^ at the cusp q = e^^'^ = (see Eq. ( |2.9| ), thus reproducing, for k = 1, the simple 
pole associated to the unphysical tachyon of the heterotic string. The set of Niebur-Poincare 
series J-'(s, k,w) is closed under the action of the derivative operators Da, and Da, defined in 
( |A.2[ ), which, according to ( 2.11| |, act by raising or low^ering the weight w by tw^o units while 



keeping s fixed. 

At the special point s = 1 — j, which for w < lies within the domain of absolute 
convergence, the Niebur-Poincare series F{s,k,w) becomes a w^eak harmonic Maass formFM 



In particular, unless w takes one of the special values listed in Table 2.1 it is in general not 
holomorphic. Although the values listed in the table essentially exhaust all the cases of 
interest in string theory, it is a remarkable fact that linear combinations of Niebur-Poincare 
series, with coefficients determined by the principal part of a weak holomorphic modular 
form cpa;/ are in fact weakly holomorphic, and reproduce 0-u, itself [34 J : 

\ 
'^w= E «'«'?""' ^ ^^" = r(7- ) ^ fl,„J'(l-f,m,a;). (3.4) 



Moreover, upon using ( 2.11[ | one can also relate weak almost holomorphic modular forms 



involving (up to) n powers of £2 — or equivalently obtained by acting up to n times with 
the derivative operator Djy on a weak holomorphic modular form — to linear combinations 
of J^(s, K, w) evaluated at the special points s = 1 — j + n', with <n' < n. 

In the cases relevant to heterotic string threshold corrections, the elliptic genus 0^ has a 
simple pole at q = 0, corresponding to the unphysical tachyon, and therefore the expansion 
( |3.4| | includes only one term, with k = m = 1 (modulo an additive constant in the case 
w = 0). Moreover, the weight w is related to the signature {d + k, d) of the Narain lattice by 
w = —k/2. Since string theory restricts the Narain lattice to be even and self-dual, so that 
^d+k,d is CO variant under the full modular group T = SL(2,Z), the possible values of w are 
ly = (corresponding to the point of unbroken Es x Es or SO(32) gauge symmetry), w = —4 



^^For a definition of weak harmonic Maass forms see Section|2.3| 
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(corresponding to the point of unbroken Eg symmetry, with arbitrary Wilson Hnes for the 
other Eg factor), orw = —8 (corresponding to generic values of the Wilson lines in Eg x Eg or 
SO(32)). The complete list of weak almost holomorphic modular forms with a simple pole 
at ^ = and modular weights w = 0, —2, —4, —6, —8, —10, together with their expressions 



as linear combinations of Niebur-Poincare series, can be found in Table 3.1 Although string 



theory applications only require k = 1, our methods apply equally well for arbitrary positive 



integer values of k, which we therefore keep general until Section 3.5 



3.2 One-loop BPS amplitudes as BPS-state sums 

Since any weak almost holomorphic modular form of negative weight can be represented 
as a linear combination of Niebur-Poincare series, for the purpose of computing integrals of 
the form ( |1.2| | it suffices to consider the basic integral 



'Id+k,d{G,B,Y;s,K;T) =Id+k,d{s,K;T) = / d/^ rrf+fc,d(G,B,y) J"(s,k, -|) , (3.5) 

w^here the modular weight w = —k/2 of the Niebur-Poincare series is determined, via mod- 
ular invariance, by the signature of the Narain lattice. In order to regulate potential infrared 
divergences, associated to massless string states, we have introduced in ( |3.5| l an infrared 
cut-off T, w^hich w^e shall eventually take to infinity. 

According to the unfolding procedure, extended in the presence of a hard cut-off T in fs], 
the truncated fundamental domain J^-j- can be extended to the truncated strip {0 < T2 < 
7", — I < Ti < 2 } at the expense of restricting the sum over images in the Niebur-Poincare 
series to the trivial coset, and subtracting the contribution of the non-trivial ones integrated 
over the complement J-" — J-7-. In equations 

/•■^ dT9 /■l/^ 
Id+kA^, K,T)= -4 dri r,+,,, M^ _. (-KT2) e-2— 1 

k tI J-1/2 ' ^ (3_^) 

- 1^_^ d^rd+k,d (j-(s,K,-|)-A^^_.(-KT2)e-'^"'^^i^ 
Using the asymptotic behaviours 



el t d+k 



M^_k{-KT2)^T2\ rd+k,d-^2 ' ' as T2^0, (3.7) 



2 

and 

d 



^d+k,d ~ ^2" as T2 ^ 00 , (3.8) 



together with the Fourier expansion ( 2.14| |, one can show that the second integral in ( |3.6 



converges for !R(s) > \{2d + k), while the first integral in \3.6) converges for !R(s) > 1 + 
\{2d + k). For !R(s) in this range, one may then remove the IR cut-off and extend, in the 
first integral, the T2 range to the full R-|-. Moreover, the Ti integral vanishes unless the lattice 
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w = 



%^ = J'(2,1,0)-5J'(1,1,0)-144 

^= I T{3, \, 0) - 4J-(2, \, 0) + 13 T{\, \, 0) + 144 

^= jf5J-(4,l,0)-|j-(3,l,0) + f J-(2,1,0)-17J-(1,1,0)-144 

^ = T^ -F(5, 1, 0) - ^ J-(4, 1, 0) + I J-(3, 1, 0) - f J-(2, 1, 0) 






+ f J-(l,l,0) + i|^ 

J52g925-^(7.1.0) - ^15-^(5,1,0) + ^J-(4,l,0) - f J-(3,l,0) 
+ fj-(2,l,0)-fj-(l,l,0)-l|4 

w = -1 



A 



^J-(3,l,-2)-ij-(2,l,-2) 



¥ = 525-^(4. 1' -2) - 20-^(3. 1' -2) + M-^(2. 1' "2) 

^ = n760-^(5. 1. -2) - 350-^(4. 1. -2) + 2io-^(3, 1, -2) - ^^^2, 1, -2) 



E2 _ 1 



3 J-(7, 1, -2) - 1^35 J-(5, 1, -2) + 5I5 J-(4, 1, -2) - g^ j-(3, 1, -2) 



19819800-' V'-^' ^y 12936 -' W' -^' ^; 1 525'' \^' ^' ^> 56- 
B-^(2,l,-2) 



15- 

W 



A ~ 2520-^'. '' / 120-^(3/1/ 4j 
^i^ ~ 70560-^ (5/ 1/ —4) — 2520-^ \^' ^' ~^) + 280-^ W/ 1/ —4) 



E^ 



1 — 148648500 -^ V ' ' ' 129360-^ W'l' 4j + gggg-r (4, 1, 4j g^Qy-(3, 1, 4) 

w = — 6 



1920-^ W,l/ 6) 10080-^14,1/ 6) 



£2^4 

A ~ 241920-' V-^'-^' ^; 10080 

A 792792000-^ V'^' "/ 887040 -' V'^' -^' "/ 1 50400 



92000-*^^'^' ") 887040-*^ W/1, 6) + gg-Qgy- (4, 1, 6j 



w = -8 


4" - 


2854051200-^('^'l' ^) 3991680-^(^' 1' 


-8) 






w = -10 




1 



1 = 131-^(7.1.-10) 

Table 3.1: List of all weak almost holomorphic modular forms of negative weight with a 
simple pole a\.q = Q, as linear combination of Niebur-Poincare series J^{\ — '^+n,\,w) (the 



holomorphic ones appear in the first column of Table 2.1 1. 
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vector satisfies the level-matching constraint 

pi-pi = AK. (3.9) 

In heterotic string vacua (with k = 1) this condition selects the contributions of the half-BPS 
states in the perturbative spectrum, and thus the first integral in ( |3.6| can be written as a 
BPS-state sum 



IduAs.K) ^ L r^A^,,_a-'<T-2)Tf2e— 2(P^+Pi)/2. (3.10) 

RPC JO Tr. 2 



dT2 

BPS -"^ '^2 

Here we have introduced the short-hand notation 



BPS PL , PR 



to denote the sum over those lattice vectors satisfying the level-matching condition ( |3.9| l, and 
corresponding to half-BPS states if k = 1. By the previous estimates, this sum is absolutely 
convergent for K(s) > 1 -|- 5 {2d + k), and thus defines an analytic function of s in this range. 
To relate the BPS-state sum to the modular integral of interest, we note that upon using 
( |3.6| l and rearranging terms, Eq. ( 3.10| maybe rewritten as 



^d+k,d {S, k) =Xd+k,d (S, K, T) 

+ l^_^ dlir,+^,, ('j-(s,K,-|)->^^_|(-KT2)e-2--i _/0(s)t2"^+^') 

+ / dll{rM,,-4) fx,_,(-KT2)e-2-'-i+/o(s)T2 



1-s+l 



(3.12) 



where 



(47r)i+l7r^i2 r(2s - 1) k^+i (Tx_2,{k) 



his) = ^ — '- r^ 1 izfi^^ (3.13) 

^^^ r(s + 1) r(s - I) a2s) 



is the coefficient of the zero-frequency Fourier mode ( |2.15 l, and the rh.s. of ( |3.12| | is inde- 



pendent of T. The first three lines in ( |3.12[ | are analytic functions of s for !R(s) > 1, since 
Xd+k,d{s, K, T) is integrated over the compact domain J^-j-, while the integrands in the second 
and third line are exponentially suppressed as T2 — > 00, away from the points of enhanced 
gauge symmetry. The fourth line, however, evaluates to 

/0(S) / dT2T2' ^+ * =MS)-'-^^ (3.14) 
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and is therefore analytic in s, except for a simple pole at s = \{2d + k). We thus conclude 
that the BPS state sum ( |3.10[ ) admits a meromorphic continuation to K(s) > 1, with a simple 
pole at s = ^^^ with residue /o(^^^^)- Moreover, taking the limit T — > oo in ( 3.12| |, w^e find 
that the BPS-state sum ( |3.10| l is actually equal to the renormalised integral 



R.N. / dfi rd+k,d^{s,K, -I) = \im 



'^d+k,d{s,K,T)+fo{s) 



^ 2d+k 



(3.15) 



'^d+kA^,K) 



for generic values oi s ^ ^^^- At the point s = ^^^, the renormalised integral is instead 
equal to the constant term in the Laurent expansion of 'Id+k,d (s, k) around s — ^^±^ 



4 ' 



R.N. J d;^rrf+fc,rf7"(s,K, 



lim 



X 



d+k,d 



1-d+k,d 
Id+k 



^.'^.r -/oP^ iogr+/o, 4 



■/ I Id+k 




,K 



(3.16) 



where /o(s) = d/o/ds, and the r.h.s. is defined as the limit of Ici+k,dis,K) after the pole is 
properly subtracted. 



X 



d+k,d 



Id + k 



-,K 



lim 

.,2d+k 
^ 4 



'^d+k,di^'^) 



/o 



2d+k 
4 



2d+k 
4 



(3.17) 



Eqs. ( |3.15| l and ( 3.16 1 relate the renormalised integral to the BPS state sum ( 3.10| l, or to its 
analytic continuation whenever !R(s) > 1. We note that this renormalisation prescription 
amounts to subtracting only the infrared divergent contribution of the massless states, un- 
like other schemes used in the literature where the full contribution of the massless states 
is subtracted. Of course, any two renormalisation schemes differ by an additive constant 
independent of the moduli. 

Having discussed the analytic properties of the BPS-state sum ( 3.10| l, and its relation to 
the regulated integral ( |3.5| l, let us now evaluate the integral in ( |3.10| |. Using the relation 
( A.16| between the Whittaker M-function and the confluent hypergeometric function ifi, as 
well as the identity 



/■OO 

Jo 



(3.18) 



we arrive at our main result 



Id+kAs,K) ={AnKy-2 r(s + ?^ - 1 



V-. „ , k 2d + k 4k 

X 2^ 2f 1 ( s - - , s + — - — - 1 ; 2s ; 

BPS 



Pi 



Pi 

4k 



2, l-s-2^ 



(3.19) 
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The sum in ( |3.19[ ) converges absolutely for K(s) > ^^^ and can be analytically continued 



to a meromorphic function on K(s) > 1 with a simple pole at s 



2rf+/c 



|34|. Again, for 



K = 1 the sum in ( |3.19| can be physically interpreted as a sum of the one-loop contributions 
of all physical BPS states satisfying the level-matching condition 10.9) . This expression is 
manifestly invariant under T-duality, independent of any choice of chamber, and generalises 
the constiained Epstein zeta series considered in 1 20 2T| to the case of a non-trivial elliptic 
genus. We would like to stress that these properties follow^ directly from our approach, as 
opposed to the conventional unfolding method, which depends on a choice of chamber to 
ensure convergence. 

Moreover, using the fact that the lattice partition function satisfies the differential equa- 
tion 1211 



A 



SO(d+k,d) 



2A^k/2 + ld{d + k-2) ra+k,d = 0, 



we find that the BPS state sum ( 3.19| is an eigenmode of the Laplacian Aso((i+)c,rf) 
Narain moduli space 



A 



For s 



SO{d,d+k) "I" 16 



2d+k 



j^(2d + /c-4s)(2d + /c + 4s-4) Xrf+^,rf(s,K) = . 



(3.20) 
on the 

(3.21) 



, the eigenvalue vanishes but the BPS state sum X^^j^ii{s,k) has a pole. After 
subtracting the pole, one finds that the renormalised BPS state sum is an almost harmonic 
function on the Narain moduli space, namely its image under the Laplacian is a constant 



Id + k 

AsOid,d+k) -l-d+k.d 1 -, ,K 



(l-rf-|)/o 



2d + k 



(3.22) 



3.3 One-loop BPS amplitudes with momentum insertions 

Our method carries over straightforw^ardly to cases w^here insertions of left-moving or right- 
moving momenta appear in the lattice sum, i.e. to modular integrals of the type 



T 



dji 



-A/2 



Y2 p{PlV^,PrV^) ^iPl^*''^ 



Pl,Pr 



cD(t), 



(3.23) 



considered for example in |26l38|. The term in the square bracket is a modular form of 
weight (A + d + 2/0)/ provided that the function p{xi^, xr) satisfies 



[9: 



2 -92 



(3.24) 



2tz (xl3xl - ^r 9xr - a - d)] p(xl , xr) = , 

upon choosing p = T2 p^ ~ ^ (respectively P = ^2 Pr — J^)/ it is proportional to the modu- 
lar derivative D • r^j+jc,^ (respectively, D • T^j^^^,}) of the usual Narain lattice partition function. 
The integrand in ( 3.23| l is then modular invariant provided A + d + | = —w. 



and that p(xl, xr) e "^i^L+%) should decay sufficiently fast at infinity 1271^ For example. 



12 



We are grateful to J. Manschot for pointing out this reference. 
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As usual, expressing the elliptic genus as a linear combination of Niebur Poincare series, 
one is left to consider integrals of the form 



J d}iT^^'^ ^ p{VL^/Tl,pR^/T2) q'^P'^ q'^^^^ J^{s,K,w) , 



(3.25) 



Pl,Pr 



and following similar steps as in the previous subsection, one finds the result 



Bpg JO 



2d+k 



Fi s 



in 



2A + 2d + k 



;2s;H p 



PL 



v/4 



VI 



Vr 



7ZK 



Vi 



TCK 



(3.26) 



In most applications, p is a polynomial in p^/ Pr/ ^rid the integral can be evaluated using 



(3.18 1. As a result, each monomial can be evaluated to 



Pl,Pr 



/ d^ 4 E Pl • • ■ Pl Pr--- Pr q'^^'q''^'^ J^{s,K,w) ^ {AnK^-^ r(s + ^ + ^ - 1) 

\-^ fli ah, ba -^ f \w\ \w\ ,- , ^ 4k:\ /p?\ 



i_s_m_^ 



(3.27) 



with (5=(a + /5 — A)/2. Clearly, this result is meaningful only when the various monomials 
are combined into a solution of ( |3.24| , as required by modular invariance. 

3.4 BPS-state sum for integer s 

For special values of s and w, the hypergeometric function 2^1 appearing in the BPS-state 
sum ( |3.19[ ) can actually be expressed in terms of elementary functions. For example, for 
d = landw = 0, 2Fi(s,s - i,2s;z) = 22«-i(l + ^l-z)^-^', and thus 



X14 (1 + n, k) = VAtzk 2^+2" r{n + l)Y^ 



BPS 



4k V 4k 



-l-2n 



lv/7r(16K)i+"r(n + l) X: 






pR + qR 



-1 



pR-qR 



-1 



-l-2n 



(3.28) 



w^ith s = 1 + n. For n = 0, this agrees w^ith the expression derived in |20 1 using the Selberg- 
Poincare series S{s,k,w) at s = 0. 

More generally, similar simplifications also take place for s = 1 — ^ + n = l + |+n, 
w^ith n a positive integer, which are the special values relevant for representing w^eak almost 
holomorphic modular forms, and are thus of interest for our physical applications. While 
it is cumbersome to express 2F1 directly in terms of elementary functions, it is simpler to 
notice that the Whittaker M-function appearing in ( |3.10| reduces to the finite sum ( |A.33[ ). As 
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a result, the integral ( 3.10| l reduces to 

RPQ >^0 4 ' 2 



BPS 



(47rK 



i_^r(2(n + l) + §)r(n + ^^ " 



n! 



m=0 



m 



(-1)' 



r{n-m + ^^ 



d+k^ 






(3.29) 

where, in going from the first to the second line we have set z = 4lTckt2, and have integrated 
by parts n times, and we used the fact that the boundary terms vanish. Although the full 
integrand vanishes rapidly enough as z — > 0, so that the integral exists, this is not true of 
each individual term, unless ^ — n — 1 >0. To regulate these unphysical divergences, we 



introduce a convergence factor T2 in the integrand, and evaluate each integral in (3.291 for 
large enough a. The desired result is then expressed as the limit 



ld+kA'^ + ^+n,K) =(47rK)^-3 



1 . r(2(n + l) + |)r(n + 



L 

BPS 



9 X n—m 

Pi 

4k 



lim 



d±k\ n 
m=0 



ri I -m-l + OL 



{-ly 



mj r{n-m + ^] 

pV 

4k 



L 



2n+k/2 r {^-m-l + i + a 



4k 



(3.30) 



The series ( 3.30| converges absolutely for n > | — 1, as a result of the finiteness of the orig- 
inal modular integral. For n < ^ — 1, it is a formal (divergent) sum over BPS states, which 
nevertheless captures the singularities of the amplitude at points of gauge symmetry en- 
hancement. 

For n < 2 — 1/ or whenever d is odd, independently of n, the limit a — ?> is trivial, leading 
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(1) 



to ld+kA^> k) = Tl^^^s, k) where 



d/2-2 

X E 

m=0 



1 .r(2(n + !) + §) r(n + ^ 



d+k\ 



n\ 



m 



{-If 



r{n-m + ^ 



E 



d+k\ L^ \ A^ 
BPS ^ ^'^ 



7 X n—m 
PI 



r[{ 



m 



4k 



2n+)c/2 r(|-m-l+^ 



E 



£! 



(3.31) 



4k 



If d is even and n > ^ — 1 one finds Id+k,d{s,K) = '^d+kdi^'^) +-^i+*:d(^''^) where the first 
term is still given by ( |3.31| l and the second term is 



j(2) n I k 



-d+kM 



(1 + 



. ,^ ,i_rf r(2(n + l) + |)r(n + ^ 
n,K) = {AtikY 2^ / 2^ V 2 



xE E 
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my r{n-m 

d 



d+k^ \^4^ 
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in+k/1 r[^-m-l + £ 

i=m+l-d/l 
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4k 
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(-ir 



r{m 



2- {) V4k 



2 X m+l-'j 
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m+l- 



.OS ( I 



m+l-rf/2 



■r(^ + 2-f) S 



m + 1 



Pl 
4k 



2 ^ m+l-^-i 



H 



m+i-^-e 



> , 

(3.32) 



where H^ = I^k=i ^^ i^ ^^^^ N-th harmonic number. The combination ( 3.31| l vanishes for 
d = 2, the sum over m being void. The results ( |3.31 1 and (3.321 allow to write any integral 
of the type ( |1.2| as a formal sum over physical BPS states (w^hich converges absolutely for 



n > 2 — 1). In particular, the result is manifestly invariant under the T-duality group 0(d + 
k,d;Z) 

We conclude this subsection with some simple examples for special values of n and k. 
For n = the sum over m in ( 3.29| l is void and only few terms contribute to the integral, 
corresponding to the various terms in ( 2.23| |. When d ^ 2 the limit a — ?► is trivial and one 
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arrives at the simple expression 



BPS 



nf-l)(g 



i—l p\ 



11=0 



4k 



(3.33) 



When d = 2, the Hmit a — > is subtler and leads to logarithmic contributions. One obtains, 
for n = 0, any k, 



i2+k,2{i + lK) = -r{2 + l)Y: 

BPS 



^<ym^ 



(3.34) 



and for k = 0, any n, 



I2,2{1 +n,K) = ^^ l^ 



„2x " n 



n! ^ V4k / '— '„ \m 

BPS \ ^ m=0 



Pi 



H,„-log{4 
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E(-i) 

^=0 



£/"^\ (Pl 
£/ Uk 






^! 



4k 



(3.35) 



As usual, the left and right-handed momenta are defined by 

Pl,i = {mi + Y^Q" + lY^Y^ni + (G + B)ijni , Q« + Yfri') , 
PR,- = WJ; + Y^Q" + lyf Y;«n^' - (G - B)ijni , 



(3.36) 



with m, and n' the Kaluza-Klein and winding numbers and Q" are the charge vectors. In the 
d = 2 case, and in the absence of Wilson lines, it is often convenient to express them in terms 
of the Kahler modulus T and of the complex structure modulus U as 



pi 



Pl 



T2U2 



T2U2 



m2 — Umi + T{n^ + Un^ 
m2 — Umi + T{n^ + UrP' 



(3.37) 



The relation between these results (for A: = 0) and the 'shifted constrained Epstein zeta 
series' of |20J is discussed in Appendix |BJ 

3.5 Singularities at points of gauge symmetry enhancement 

In addition to keeping T-duality manifest, another advantage of this approach for the evalua- 
tion of one-loop modular integrals is that it allows to easily read-off the singularity structure 
of the amplitudes at point of enhanced gauge symmetry. These points are characterised by 
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the appearance of extra massless states with pR = 0. Depending on the dimension of the 
Narain lattice, as well as on the value of n, the amplitude may diverge ( we refer to this case 
as real singularity) or one of its derivatives can be discontinuous (w^e refer to this case as 
conical singularity). 

For odd dimension d, the modular integral Id+k,d{S' t<^) always develops conical singular- 
ities, as exemplified in the one-dimensional case by Eq. ( |3.28| . In addition, for d > 3 real 
singularities appear from terms with m < | — 1 in ( |3.30[ ). 

For even dimension real singularities alw^ays appear. They are are power-like in Z^^> 
w^henever d > 4 and logarithmic in X^^ for any even d < 2n + 2. Moreover, conical singu- 
larities do not appear. 

Notice that for d = 2 the singularities cancel in the combination 

X2,2(l + n,l)- ^ ^; ' X2,2{l,l), (3.38) 

which is therefore a continuous function over the Narain moduli space, including at points 
of enhanced gauge symmetry. Since, using the results in | 10p8 | and the fact that F{1, 1, 0) = 
7 + 24, 

^2,2(1,1) = -loglK^) -7(^)l'-241og \T2U2\r]{T)ri{U)\^'\ +const, (3.39) 

we conclude that all integrals 12,2(1 -\-nA) exhibit the same universal singular behaviour up 
to an overall normalisation 

X2,2(l + n,l) ~ -^^^^ log|;(T) -7(^)1'- (3.40) 



This expression can be generalised as in 1 18 1 if Wilson lines are turned on 



4 Some examples from string threshold computations 

In this section w^e evaluate a sample of modular integrals that enter in threshold correc- 
tions to gauge and gravitational couplings in heterotic string vacua using the method de- 
veloped in the previous section. We express the elliptic genus as a linear combination of 
Niebur-Poincare series, and w^e evaluate the modular integral in terms of the BPS-state sums 
^d+k,d{s,K) defined in eqs. ( 3.19 1 and ( 3.30| |. 



4.1 A gravitational coupling in maximally supersymmetric heterotic vacua 

Let us start with the example of toroidally compactified SO(32) heterotic string, for which 



the elliptic genus takes the form ( |3.1| |. Using Table 3.1 and the relation E^A ^ = 7 + 744, this 
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(4.1) 



can be conveniently expressed in terms of the Niebur-Poincare series as 

0{t) =ts tr F4 + ^^ [ J-(l, 1, 0) + 720] h tr R^ 

+ ^ [- i J-(3, 1, 0) + 5 J-(2, 1, 0) - 18 J-(l, 1, 0) + 288] ts tr R^ tr F^ 
+ ^ilJ'iXhO) -6J^{2,1,0) +2AJ^{1,1,0) -576] tsitvF^f. 

Therefore, using the results in the previous section, the renormalised modular integral ( |1.2| | 
can be expressed as the linear combination 

R.N.| d^^d40=Id,dtstrF^ + ^^ [Xrf,rf(l,l) + 720Xrf,rf] fgtrR^ 

+ ^[-lldA^A)+5IdA^A)-18IdA^A)+288Id,d] tstvRhrF^ 



^9^ [iXrf,d(3,l)-6Xrf,rf(2,l) +24X^^,^(1,1)- 576 Xd,d] tsitrF^f. 



(4.2) 



where, as computed in 1 20 21], 



Xd,rf = R.N. y^d/^r^,d(G,B) 



r( 



^y(G, B; 2 — 1), 



(4.3) 



7l2' 



with £y being the constrained Epstein zeta function defined in 1 20 21 . In this expression, 
any time n = j — 1 the BPS-state sum X^A^ + ^' '^) should be replaced by Xd,i(l + n, k), as 



explained in Section 3.2 



In the one-dimensional case the constrained sums can be easily evaluated, leading to 

I d^ri,ia)=|(R + R-i)f8trF^ + ^3^(l5R + 16R-i) fgtrR^ 



71 



2532 

n 
¥3 



3R + 16R"^ - 24R-3 + 12R-5 ) fg (trR 



72\2 



R - 2R"^ + 5R-^ -2R-^) f g tr R^ tr F^ 



(4.4) 



-^{R-R-'+3R-'-R-')tsitrF')\ 

for R > 1. The expression for R < 1 can be obtained by replacing in the previous expression 
R I— > R^^. Notice that, aside from the threshold correction to tg tr F'^, all other terms develop 
a conical singularity at the self-dual radius R = 1. 
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4.2 Gauge-thresholds in J\f = 2 heterotic vacua with/without Wilson lines 

Let us turn now to A/" = 2 heterotic vacua in the orbifold Hmit T^ x T^/'Z.2, with a standard 
embedding on the gauge sector. At the orbifold point, the gauge group is broken to 



Eg X Eg ^ Eg X E/ X SU(2) , 
and, in the absence of Wilson lines, gauge threshold corrections read 

^E« = -^y^dHr,,. K 

W E2E^Ee-El 



^E. 



From Table |3J] one can read that 

J'(2,l,0)-6J'(l,l,0)+864, 



El £4 Es — £5 



A 



E2 £4 E(, — E4 



J'(2,l,0)-6J'(l,l,0)-864, 



(4.5) 



(4.6) 



(4.7) 



and thus 
A 



L 

BPS 



E7 - I] 
BPS 






+ 72 log (T2U2 \i]{T)rj{U)n + const. 



72 log [T2U2 \i^{T)}]{U)n + const 



(4.8) 



Notice that the combination ^2(2, 1,0) — 6X2(1, 1,0) is regular at any point in moduli space 
(and in any chamber), as expected since the unphysical tachyon is neutral and therefore does 
not contribute to the running of the non-Abelian gauge couplings. 

Turning on Wilson lines on the Eg group factor along the spectator T'^, yields 



1 /■ , r ^^Ee-El 



(4.9) 



Using Table 3.1 one easily finds 



^"^^ = It{A, 1, -4) - ^ J-(3, 1, -4) , 



(4.10) 



and thus 



Ae, 



1 

720 



42 



Xio,2(4,l)-Xio,2(3,l) 



E 

BPS 



16 64 



1 + PiiogfPiV--^ 



(4.11) 
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In this expression pl,r depend also on the Wilson lines, and the constraint in the BPS-sum 
now reads 

PL-Pi = 4 => m'n + lQ'Q = l, (4.12) 

where Q is the U(l)-charge vector in the Cartan sub-algebra of Eg. 

4.3 Kahler metric corrections inAf = 2 heterotic vacua 

Our procedure can also be used to compute loop corrections to Kahler metric and other 
terms in the low-energy effective action. For instance, in A/" = 2 heterotic vacua at the 
orbifold point, the one-loop correction to the Kahler metric for the T modulus reads 



ixyy 



2 (4.13) 



-^ / d^(J-(2,l,-2)3^r2,2, 



72,. X2 



where we have used the relation between E4E^A^^ and T{s, k, w) from Table 2.1 Integrating 
by parts, and using the action of the modular derivative on the Niebur-Poincare series, one 
immediately finds 



K 



TT 



1— loop 



Similar results can be obtained for higher-derivative couplings in J\f = A vacua. 

4.4 An example from non-compact heterotic vacua 

In some heterotic constructions on ALE spaces and in the presence of background NS5 



branes, gauge threshold corrections include a contribution of the (finite) integral 1 45 1 



L = J^d,WT,rjfjr^'^^^^'^^^. (4.15) 

Despite its apparent complexity, this integral can be easily computed using our techniques. 
In fact, the relation 

!?i(0|T) =27^?/^ (4.16) 
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and the standard bosonisation formulae, allow one to write 



iV^vvf 



1 3 
'SndR 



-(ri4(2R)-ri4(R)) 



R=l/y2 



(4.17) 



Combining this observation with Table 3.1 Eq. (328), and with the standard result J-p dfi Ti i (R) 
^ (R + R^^), the integral reduces to 



1 3 
"S/r 3R 



l^dH(ri,i(2R)-ri,i(R)) 



X (^-F(3, 1, 0) - 6J'(2, 1,0) + 23 J'(l, 1, 0) + 432) 



R=l/V2 



1 3 
"Stt 3R 



^ ^ 1 Xi4 (2R; 3, 1 ) - 6 X14 (2R; 2, 1 ) + 23 X14 (2R; 1, 1 ) 



;iXi,i(R;3,l)-6Xi,i(R;2,l)+23Xi,i(R;l,l))+144 7r(R-lR-i) 



-20 v/2 . 



(4.18) 
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A Notations and useful identities 

In this appendix we collect various definitions, notations and formulae used in the text. 



A.l Operators acting on modular forms 

The hyperbolic Laplacian acts on modular forms of weight w vi-J^ 



A. 



iT^df [dr 



2^j 



(A.1) 



^^Our Laplacian is related to the one used e.g in 36 via A^ = — 2 An);BO 
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We denote by 'H{s,w) the eigenspace of Aiy with eigenvalue 55(5 — 1) — \w{w + 1), in the 
space of real analytic functions of modular weight w under F = SL(2, Z). The raising and 
lowering operators D^, D^ defined by 

D^ = ^ {d^ " ^) ' ^"' = "'"^ ^2 9t , (A.2) 

map 'H{s,w) to H{s,w ±2), 

n{s,w-2)^n{s,w)^n{s,w + 2), (a.3) 

and satisfy the commutation identity 

- - w 

D,„_2 ■ D,„ - D-a,+2 • O™ = - . (A.4) 

The operator D^, (and of course. Da,) satisfies the Leibniz rule 

Div+zo' [fiu fw') = {Dw f-iv) fiv' + fm (Da,' /a,') , (A.5) 

where /a, is a modular form of weight w. We denote by DJ„ /a, (or simply D*"/) the iterated 
derivative Da,+2r-2 ■ ■ ■ ■ • Da,+2 • D^, ■ f^,, a modular form of weight w + 2r. One has 

For a; < 0, the operator D},^^" simplifies to {i/7iy^^d]r^" (Bol's identity), and is known in 
the physics literature as the Farey transform p2| . 
The Hecke operators Tk are defined by 

a,d>0,ad=K b mod rf ^ ^ 

and satisfy the commutative algebra 

T,T,,= X] d'-'"T,,,/,2, (A.8) 

rf|(K,K') 

li (p = Y^„,=z ^{^f T2) e^^"^i is a modular form of weight w, then the Fourier coefficients of 
T^ ■ are 

{T„-0){n,T2) = K'-'" Y. d'"-'0{nK/d^,dh2/K) . (A.9) 

(?|(n,K) 

The generators of the ring of holomorphic modular forms are the normalised Eisenstein 
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series 



E^ = l+2A0j2r=^ ^^ £6 = 1-504X]7:^, (A.10) 

with modular weight 4 and 6, respectively. The discriminant function is the weight 12 cusp 
form 

oo 

A = qYl{l- q"f' = i^(£| - El) . (A.ll) 

The generators of the ring of weak holomorphic modular forms are E4,E(, and 1/A. The 
modular y'-invariant is the unique weak holomorphic modular form of w^eight zero w^ith 

j = l/q + 0{q), 

£3 £2 

/■ = ^ - 744 = -^ + 984 . (A.12) 

^ A A 

The ring of weak almost holomorphic modular forms is obtained by adding to £4, E(„l/A 
the almost holomorphic Eisenstein series 

£2 = £2 - — = 1 - 24 ^ Jll^ _ A. . (A.13) 

Under the raising operator D^g one has 

DE2 = \{Ei-El), D£4 = §(£6-£2£4), DE(, = El-E^Ef,, D{\/ A) =1E2/ A, {A.IA) 

where, for simplicity we have left implicit the specification of the weight in D. Using the 
Leibniz rule ( |A.5| |, this allows to compute the action of D on any weak almost holomorphic 
modular form. 

Finally the operator T^ maps the weak holomorphic modular form 'P = \/q -\- 0{q) to 
T^<P = l/q^ + 0{q). 

A.2 Whittaker and hypergeometric functions 

Whittaker functions and hypergeometric functions, more in general, are central in the anal- 
ysis of the Niebur-Poincare series and the evaluation of one-loop modular integrals. We 
summarise here their definitions and some of their main properties. 

Whittaker functions are solutions of the second-order differential equation 



""+ I -7 + - + ^^ I " = 0- (A-15) 



4 2 2^ 
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For 2}i not integer the two independent solutions are given by the Whittaker M-f unctions 

M;,,±^{z) = e-^/2 2±f'+5 ^f^ {±^-A+l;l± 2ii;z) , (A.16) 

and are expressed in terms of the confluent hypergeometric function 

r(b) ^ r(a + n) z" 
r(fl) ^^rib + n) n\ 

When 2}i is an integer, however, the second solution is not defined anymore, and thus it is 
useful to introduce a second Whittaker function defined by 

W^,^{z) = -^T{X + \-^)e--'^z^ f^"^\-t)-^-l+l^ (l + i) " \-Ut, (A.18) 

•^ TTl J CO \ Z J 

w^here |arg(— f ) | < n, and the contour does not contain the point t = —z and circles the ori- 
gin counter-clockwise. The functions Wx^^{z) and W_a,^(— z) are then the two independent 
solutions of the differential equation ( A.15[ ), and have the asymptotic behaviour 



Ws,-a,(y) ~ |47ry|f(^8''(^)-i)e-27TM as \y\^oo. (A.19) 

The Whittaker M-function can then be expressed as the linear combination 

Ma„(z) = J^^^ + ^1 g'"^ W_^,,(.-z) + J^^f'^^^ g^"(^-^'-^) W,,,(z) . (A.20) 

Using the symmetry of the W-functions, Wx^fi{z) = W;^_fi{z), one can invert the previous 
relation and write 

W,,,(z) = .[^~^^\. Ma,,(z) + ./_^^^^ ^. M,,_,(z) . (A.21) 

This implies that the functions A4s,w and yVs,iv obey 

For special values of A and fi, the Whittaker functions reduce to elementary functions or 
to other special functions. This derives from the properties of the hypergeometric functions, 
for instance 

iFi(fl;fl;z) =e\ 

ifi(l,fl,z) = (fl-l)zi-''e^7(fl-l,z), 
iFi(fl,fl + l,z) = a{-z)-' 7(fl, -z) , (^-^^^ 

ifi(fl;2fl;z)=e^/2 (lz)^"''r(fl + l)J„_i(z/2), 
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where 

7(fl,z) = / e-U''-^dt = r(a)-r(a,z) (A.24) 

Jo 

is the incomplete Gamma function, and 

is the modified Bessel function and /i/(z) is the Bessel function of the first kind. 
Given the definitions ( |A.16| |, ( |2.7| | and 



Ky{z) 



n I-v{z) -Iv{z) 



2 sin TTV 
one can then show that (y > 0) 

Ws,o{±y)=2\y\'^K.{2n\y\), 



(A.26) 



(A.27) 



>(u,„r-y) = e2-y. 



W|^^„(-y)=^(l-^(;,47^y)e2-^ 



Wi. (y) = e 



-Iny 



(A.28) 



M 



i-y) = {^nyre~^-y, 



W_^i-y) = {4.nyr"e 



w p-lny 



y^^^^Sy) = {^ny)-^T{\^wAny)^^\ 



W^_-^,„{-y) = r{l-wAny)e'^y , 



Wi_|,,„(y)=e-2-y. 



(A.29) 



(A.30) 



For integer values of the arguments, the confluent hypergeometric function, and thus the 
Whittaker functions, take a particularly simple expression 



ifi(n + l,fl,z) 



^n+l-a . g 



r(fl) d" 

n\ dz" 

r{a)z'-^\e^Lt"\-z)-Lttni-) 



a-l ^k 

z Y^ _ 

k=0 '^- , 



(A.31) 
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where 



L^W 



n! dx" 



^n+k^-x 



(A.32) 



are the associated Laguerre polynomials. 

As a result, the seed function that enters in the definition of the Niebur-Poincare series 
involves only a finite number of terms when s = 1 — |^ + n, and reads 



\—-2+n,w 



{-y) =(47ry)i-'''+"e-27ry 
d" 



(2n + l -w)\ 



;! 



[AnyT-'^-' U-y -Y^^ 



d(47ri/)" 
=T{2n + 2-w){4:ny)-" 

Similarly, 

>^i-f +.,Jy) =(-1)""! (47ry)-".-2-^Ll-i-^"+^'')(47ry), 



(A.33) 



(A.34) 



The modular derivatives ( A.2| | have a natural action on the Whittaker functions, so that 






(A.35) 



and 



D^ ■ [W.,^(nT2) e^-"-!] = W.,,.+2(nT2) e^'^^''^^ 



D.- [W.,.„(nT2)e2"^"^i] = W.,.„-2(nT2)e2--i 



-2n, 



n >0, 



2n(s + f)(s-f-l), n<0. 



8)1 



n < 0, 

l(s-f)(s + f-l) n>0. 



(A.36) 



8n 



A.3 Kloosterman-Selberg zeta function 

The Kloosterman-Selberg zeta function entering in the expression ( 2.16| of the Fourier coef- 



ficients of the Niebur-Poincare series is defined as ||43 1 
1 



Zs{a,h) 



2^\ah\ 



L. ^ X 



if ab > 0, 



oO 



hs-i{^V^^] if ab<0, 



(A.37) 
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where Is{x) and }s{x) are the Bessel I and / functions, and S{a,b;c) are the classical Kloost- 
erman sums for the modular group T = SL(2, Z), 



S{a,b;c) = ^ exp 

dG(Z/cZ)* 



27zi 



(ad + bd-^) 



(A.38) 



Here a, b and c are integers, and d ^ is the inverse of d mod c. S{a, b; c) is clearly symmetric 
under the exchange of a and b. Less evidently, it satisfies the Selberg identity 

S{a,b)c)= ^ dS{ab/d^,l)c/d). (A.39) 

d\gcd{a,b,c) 

In the special case a ^ 0,b = 0, the Kloosterman sum reduces to the Ramanujan sum 



S(fl,0;c) = S(0,fl;c) = ^ exp ( ad 

de(z/cZ)* ^ "^ 



^ d]i{c/d), 

rf|gcd(c,fl) 



(A.40) 



with f/(n) the Mobius function. For a = b = 0, S{a,b; c) reduces instead to the Euler totient 
function (p{c), and one can verify that 



L 



S(0,0;c) a^s-l] 



r.lS 



^ S(0,±k;c) _ cti-2s{k) f^^n\ 
2-1 'tis - zi->.\ K^T^I' 



with (Txin) the divisor function. Under complex conjugation, Z^ia, b) transforms as 



(A.41) 



Zs{a,b) = Zs{-a,-b). (A.42) 

The Kloosterman-Selberg zeta function defined in ( A.37| | is related to the zeta function 

(A.43) 



, , ^ S{a,b;c) 
oO *- 



originally considered in pT| and used in [[20| via 



(A.44) 



~ ( —ATT^nh^'" 

Zs{a,b) = K{^n^\ab\r-^ Y. \r(o 1 , Z{a,b;s + m) . 

^Qm\T{2s + m) 

B Selberg-Poincare series vs. Niebur-Poincare series 



In this section, w^e briefly discuss the relation between the Niebur-Poincare series ( |2.8| and 
the Selberg-Poincare series ( |2.5| l, considered in our previous work pO| in the special case 
w = 0, as well as the relation betw^een the BPS-state sum ( 3.19| l and the "shifted constrained 
Epstein zeta series" considered in [ [20| . 

Comparing the differential equations ( 2.10| and ( |2.6| , it is easily seen that a set of solu- 
tions of one can be converted into a set of solutions of the other by considering the linear 
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combinations ||33l|44 



F{s,K,w) = Y2 a{s,K,w,m) £(s + m,K,zv) , 



"2,° (B.l) 

E{s,K,w) = 2_j h[s, K, w, m) F{s + m,K,w) , 

m>0 



such that the coefficients satisfy the recursion relations 

a{s,K,w,m + l) AnK{s + m — j) b{s,K,w,m + l) 47rK(s-f) 



(B.2) 



a{s,K,w,m) (m + l)(m + 2s)' b{s + l,K,w,m) (m + l)(m + 2s) 

Comparing also the constant term l\2.15) of the Niebur-Poincare series and the constant term 



~ ~ 2^(1-'^) TTi-^" (tIk)"' r(2s + m - 1) ai-2s-2m{K) 1-s-m-f 

of the Selberg-Poincare series, we find that the coefficients are given by 

a{s,K,w,m) — (—1)'' 
h{s,K,w,m) 



(B.3) 



,„ 22«+2'«-'" (ttk)'- 2 +"' r(2s) r(s + m - f ; 



m!r(2s + m)r(s- f) 

T"-^' {tikY'^ 2 r(2s + m - 1) r(s + m 
m!r(2s + 2m-l)r(s- f) 



(B.4) 



In particular, in the limit s — )■ ^ where the summand of the Selberg-Poincare series ( |2.5| l 
becomes holomorphic, one finds £(s, k, w) = J^(f , k, w) for w >2, but 

E(f,K,a;) = J'(|,K,a;)+ £] &'„ReSs=|+m J'(s,K,a;) + ^ h,„F{^_+m,K,w) {^.S) 

for 1(7 < 0, where &„, = lims^w b{s,K,w,m) and &^ = limg^w j^b{s,K,w,m). In writing ( |B.5| |, 
we have assumed that the singularities of T{s,k,w) on the real s-axis can be read off from 
the constant term ( 2.15| l, namely that ^-"(5, k, w) is regular at integer values of s provided 



s>Oors<— |zf|/2, and has simple poles for integer values of s such that — ^ < s < 0. In 
particular, £(|',K,zf) receives a contribution (for m = \ — w) proportional to the harmonic 
Maass form T{\ — j,k,w), which is the main object of interest in the present work, but is 
contaminated by other Niebur-Poincare series lying outside the convergence domain K(s) > 
1. For example, for w = 0, v/e find 

£(0,K,0) = J'(0,K,0) + ij'(l,K,0), (B.6) 

consistently w^ith the identifications 

E{0,K,0) = T,i + 12a{K), T{l,K,0) = T,i + 24a{K), T{0,k,0) = JT^i . (B.7) 
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Moreover, the relation (B.ll between the Niebur-Poincare and Selberg-Poincare series im- 
pHes a similar relation between the BPS-state sum ( 3.19| l and the "shifted constrained Epstein 
zeta series" 



4(G,B,y;s,K) ^2^ j:{pl + pl-AKr = Up 



2\-s 



(B.8) 



BPS 



BPS 



generalising the constructio in (20\ to the case k ^ 0. Namely, using the same techniques as 
in our previous paper, one may show that ( B.8| l arises from the modular integral 



lim / df/r<;+^,rf(G,B,y)E(s,K,-|; 



r{s + ^-i 



n^ 4 



^4(G,B,y;s + ?^-l,K). (B.9) 



Since the BPS-state sum Id+k,d{s, k) arises from the limit T — > oo of the integral ( |3.15| l, from 
( B.l| | we conclude that 



TT' 






£:^(G,B,Y;s + g^-l,K)= /^ ^^^^ J2 Hs,Kzv,m)IduAs + m,K) , (B.IO) 

^ \S-\ 4 1) „,>o 

for large 3?(s). 
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